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Abstract. The Fermi two-atom problem illustrates an apparent causality 
violation in Quantum Field Theory which has to do with the nature of the built 
in correlations in the vacuum. It has been a constant subject of theoretical debate 
and discussions during the last few decades. Nevertheless, although the issues at 
hand could in principle be tested experimentally, the smallness of such apparent 
violations of causality in Quantum Electrodynamics prevented the observation 
of the predicted effect. In the present paper we show that the problem can 
be simulated within the framework of discrete systems that can be manifested, 
for instance, by trapped atoms in optical lattices or trapped ions. Unlike the 
original continuum case, the causal structure is no longer sharp. Nevertheless, 
as we show, it is possible to distinguish between "trivial" effects due to "direct" 
causality violations, and the effects associated with Fermi's problem, even in such 
discrete settings. The ability to control externally the strength of the atom- 
field interactions, enables us also to study both the original Fermi problem with 
"bare atoms", as well as correction in the scenario that involves "dressed" atoms. 
Finally, we show that in principle, the Fermi effect can be detected using trapped 
ions. 
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1. Introduction 

Fermi [1] considered two atoms A, B, coupled to an electric field. Consider The initial 
state It^is Of) of the system , where fyi and Ib denote atoms A and B in the excited 
and ground state respectivly, and 0^ - the ground state of the electric field. What 
is, then, the probability of transition to the state liAts Of)? Causal considerations 
lead us to believe that the probability must satisfy P cc d (t — This is the result 
obtained by Fermi, perturbatively, by making a further assumption, extending the 
domain of a frequency integration to include negative frequencies as well. Shirokov, 
however, showed that the result is noncausal, and the causal probability was found by 
Fermi and others due to this approximation (see, for example, his review [2]; Other 
method of calculation can be seen in O |4]). The reason for that behaviour was 
examined by several authors (such as O [Sj El [?])■ A common solution to the problem 
is to reformulate the problem, stating that such a probability is not measurable, since 
it requires a nonlocal instantaneous measurement of both the atoms (and the field). 
Therefore some authors did not post select the state of A or of ^ and the field[||. It was 
shown in various ways, either perturbatively or exactly, that the transition probability 
when only B is measured is indeed causal, in the sense that it is independent of atom 
^ for < < R/c (where R is the distance between the atoms), although it was shown 
that noncausal correlations develop between the ions (see, for example, [6J). 

This leads us to the issue of dressing: It was shown mathematically by Hegerfeldt 
[SI that within a renormalized theory, i.e. considering the dressed states and 
Hamiltonian rather than the bare ones, the probability to find atom B excited (without 
post selecting A and the field) is not proportional to 6 (t — -^). 

It is therefore reasonable to believe that the existence of a state's "tail", or of 
a "photon cloud" is the cause for what can be interpreted as a noncausal result. 
Renormalization was introduced by Hegerfeldt as a possible "way out"; The dressing 
process of a two- level atom coupled to a field was investigated [9] : It was shown that 
both in bare and dressed initial states of a single atom coupled to the field, there is a 
nonlocal photon cloud surrounding the atom. In the bare case it is canceled for t = 0, 
and hence this nonlocal behaviour is hidden in the initial condition ("curtain effect"), 
however, for any positive time there is a photon distribution around the photon. 

The Fermi Problem has been a constant subject of theoretical debate and 
discussions. Although the issues at hand could in principle be tested experimentally, 
the smallness of the apparent violations of causality in Quantum Electrodynamics 
prevented the observation of the basic effect as well as the study of the additional 
modifications due to the dressing of the photon cloud. 

The first objective of this article is to show that the Fermi problem can in fact be 
reformulated in the framework of discrete systems. This has the clear advantage as it 
enables to simulate the problem within the context of physical atomic systems such 
as trapped atoms in optical lattices [lOJ or trapped ions [TTl I12|. In another recent 
proposal, the continuous eflect has been studied in the context of circuit QED [13J. 

However, the discreteness of our systems raises a new problem: the concepts of 
"causality" and a "light-cone" now become ill defined and hence, non causal effect 
might obscure the effect we are looking for. We shall have a close look at such errors 
and show how to distinguish between dynamical non-locality due to the finiteness of 
the system and the Fermi Problem. 

f Some Authors had even used a more general initial condition on A, or both A and B. 
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The use of atomic discrete systems has however many important advantages; 
First, unlike QED, the use of vibrational degrees of freedom (phonons) as simulating 
the mediating field particles, and that of the internal electronic levels as the "atom" 
detectors, allows for different time scales for propagation and measurement, i.e. the 
measurement can be done at a much shorter time scale compared with the propagation 
time between the atoms. Second, thanks to the developed quantum information 
techniques, the interaction between the "atoms" and the "field" can be accurately 
controlled and hence turned on and off [12]. We will exploit this ability in order 
to control the dressing or the phonon cloud, or sometime to avoid the dressing and 
access bare states; something that in QED cannot be done. Finally, using quantum 
information techniques, it is possible to measure both the internal degrees of freedom 
of the "atoms" A and B, as well as the motional phonon state of the discrete field. 
Hence after turning the field-atom interaction off, before the "causal time", one 
can perform, effectively, an instantaneous, nonlocal measurement, by measuring the 
systems separately, which is now possible since they are independent. Therefore one 
can also simulate the original model envisioned by Fermi, which assumes bare atoms 
and a final post selections of the atom A and the field. 

The article proceeds as follows: In section 2 we discuss the effect in discrete 
systems. After stating the problem, we discuss the concepts of "causality" and an 
approximate "lightcone" in such systems, which are different, of course, than in the 
continuous case. We then present our study of the dressing process and discuss the 
possible scenarios in the presence of dressing. 

In section 3 we show examples of the effect in two discrete systems: the linear 
harmonic chain, and ions in a linear trap. 



2. The Effect in Discrete Systems 

2.1. Statement of the Problem 

Consider a set of N vibrating external degrees of freedom - atoms, ions etc. We 
assume that the system has a quadratic Hamiltonian. The eigenmodes are created 

and annihilated by aj,, a^, ( a^, aj = Ski), and the eigenmodes' frequencies are {cok}- 

Hence, the displacement (from equilibrium) operators have the general form 



k=0 



(A„fcafc + Xlkai) (1) 

k=0 

and the Hamiltonian is = 1) 



N-l 



Hp = ^Wfea^Ofc (2) 

fe=0 

Besides the external degrees of freedom, we assume two-level internal degree of 
freedom on each site, represented by a spin-1/2 system, with the free Hamiltonian 

i^O = \^n< (3) 

which interacts with the external (vibrational) degree of freedom in the same site in 
the following local manner: 
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Hf = ef,,it)q,,a^. (4) 

The Fermi Problem and other related problems where atoms are coupled to a 
field were mostly formulated with the atoms coupled to an electric field E, but also 
to a scalar field, where the coupling is either to its conjugate momentum (such as 
in [9]) or to the field operator itself (such as in |7]). The 'noncausal' effect should 
not differ across these cases, since field propagators are not causal, either in scalar or 
electromagnetic field theory. Therefore, for the sake of simplicity, we shall work in the 
scalar case, where the atoms are coupled to the field. Moreover, this is the interaction 
used in the case of ion traps, as in the relevant section in this paper, and hence it is 
used throughout other parts of the paper as well. 

The total Hamiltonian is then: 

H = Hp+J2{HS + Hf) (5) 

n=0 

Note that this Hamiltonian can be implemented in various ways, such as using 
ion traps (|11|, I12|) and optical lattices ([10l[T4]). One of the main requirements from 
the system is the possibility to measure the internal degree of freedom fast enough. 

After defining the physical systems of interest, we can state Fermi's Problem. 
The external motional degrees of freedom will serve as a discrete analog to the field. 
We pick two sites A and B to act as the atoms which carry internal levels. Therefore 
we set, for every n ^ {A, B}, fn (t) = 0. 

The system is prepared in the initial state 

\^{t = 0)) = \l^,) = nAiBOF) (6) 

and the transition probability P (t) to the final state 

Hf) = \UtB Op) (7) 

is examined. 

As in Fermi's original formulation, after turning the interactions off, the state of 
all three systems {A, B, F) can be measured separately, in an effectively instantaneous- 
nonlocal measurement, which is not possible if the interactions are constantly open. 
Therefore this post-selection, which was "forbidden", is possible, since it can be 
performed after the interactions are switched off. 

2. 2. Causality 

We shall develop here a method for distinguishing the a-causal contribution that has a 
dynamical source, from the a-causal contributions that arise in the Fermi Problem. Let 
us calculate the transition amplitude using the interaction-picture, with the interaction 
Hamiltonians (for A, B) 

= e/„ (t) «e'""* + a!!e-'""*) ^ [XukakC-''^"' + X^A"'""') (») 

k 
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The I'^j) — > \tpf) transition amplitude is (see appendix 1 for details) can be 
written (in leading order) as 

-^ = \ f dt'fA it') e-"^-*' I'df'fB it") e^''-'" (0| {qa {t') , Qb {t")} |0) + 
^ Jo Jo 

t t' 

\ ^ rfi'/A (iO e-'"-*' (^ d£'!B it") e'"«*" (0| [qa it') , QB (*")] |0) - 

dt"fB it") e^"-*" (0| [qa it') , QB {t")] |0)) (9) 

Hence it decomposes to two parts: the first is commutative, correlation-dependent 
and contains two independent time integrations, 



Ao (t) 1 



= \J^ dt'fA it') e-*"-*' dt"fB {t") e^f^«*" 



{qA{t'),qB{t")}\Q) (10) 



and the second is non-commutative, commutator-dependent and integration order 

dependent, 

Ac jt) ^ 

\ £ dt'fA it') e-^"-*' ij^ dt"fB it") e^''-'" (0| [qa {t') , QB {t")] |0) - 

£ dt"fB it") e'"-*" (0| [qa it') , QB it")] |0)) (11) 

In the continuum field limit, Qnit) — > 0(a;„,f), the commutator vanishes 
for a spacelike separation (t < lEs^JEAl'j and therefore only the contribution of 
Aq it) remains. This is the part which dominates and "causes" the nonlocal Fermi 
effect; Therefore, we would like to show that for discrete systems, where Ac it) is 
nonvanishing, its contribution is much smaller or negligible compared with Aq it), for 
times that are restricted within the effective lightcone. 

Generally, 

Qn it) = ^ (A„fcafee— --^ + A;,4e-''*) (12) 
fe 

and hence, defining t = t" — t' , 

(0| QA it') QB it") |0) = J2 (0| AAfeA^feafe4 |0) = ^AAfeA^fee-'=- (13) 
fe fe 

Define the functions 

Fa (r) ^ (0| {qA (to , QB in} |0) = ^ (AAfeA^^e-''- + A^feABfee--^'^ (14) 

fe 

and 

iFc (r) = (0| [qA it') , qB it")] |0) = ^ (AAfeA*B,e^-'=- - A^^ABfee"-"-) (15) 
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In the continuum, Fc (r) ^ 9 yr — l£^_£AL j ^ So the goal is to show that in the 
relevant discrete systems, there is a typical velocity c, for which Fc (r) is negligible , 
and is particularly small compared to Fa (r), for times satisfying t < liS^jSAi 

2.3. Dressing 

The diflFerence between bare and dressed states plays a major role in the Fermi 
problem, essentially because when the interaction between the field and the atoms 
is constantly open, one should not consider the bare states, but rather the "true" 
states which are the dressed ones. Hegerfeldt, in his theorem [8j, referred to the states 
as "renormalized" ones, and suggested the renormalization - the presence of a photon 
cloud around each atom - as the cause of what appears as a noncausal result. We 
interpret the "renormalization" - i.e., the process of obtaining the "physical" states of 
the system - as dressing (unlike the meaning of renormalization in QFT) and perform 
it perturbatively, as in [9l [15] . 

In the cases discussed here, however, dressing is not mandatory. As remarked 
earlier in this paper, it depends on the opening manner of the interaction functions 
/„ (t) of A and B. 

Another related issue is the field excitations distribution, which is discussed in 
appendix 2. 

Consider two different opening schemes: 

(i) The opening functions satisfy Ja.b (t) 7^ only for a very short period of time, T. 
If we set T to be smaller than the time of propagation between the two sites, we 
can be sure that positive probability is due to the effect. We can also choose T to 
be small enough so that dressing effects become irrelevant. The physical process 
can then be described with bare states. 

(ii) Turn the interactions on adiabatically, causing the system to dress. Then, at 
t = 0, bringing them to the desired dressed initial state, and allowing the system 
to evolve until t — T, when the interactions are set off. As in the previous scheme, 
T can be chosen to be smaller than the causal time. 

After turning the interactions off we have here too the benefit of effectively carrying 
out a nonlocal measurement. Let us discuss the second opening scheme in detail. 

2.3.1. Ground State Dressing Suppose we use interaction functions of the form 



Assuming, as before, that e is small enough, we can calculate the time evolution in 
time-dependent perturbation theory. We take r to be large enough, so the interaction 
is turned on adiabatically, and thus according to the adiabatic theorem. 




(16) 



we start, at t 



00 with the initial bare ground state 




(17) 



IV (t - 0)) = |G) = IgW) + e |g(i)) + Ig^^)^ + O (e^) 



(18) 
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where the corrections [G'"^) are calculated using time- independent perturbation 
theory. This is the dressed ground state. 
The first order correction M is 



(19) 



As a first order expression, it results from one operation of the interaction Hamiltonian, 
and therefore contains (as can be seen) only separate contributions from interactions 
with either A or B. Therefore, in the first order, the correction to the ground state of 
two sites coupled to the field is the same as coupling each of them separately. 
The second order correction can be decomposed into two parts: 



G^^ 



where 



E 



Bk^Bl 



kd^k 



^ V2{n + UJk)uJk 

is due to separate dressing processes {A, B alone), and 



{fl + UJk) (Wfc + Ul) 



(20) 



(21) 



^^*Ak^Bk + >^*Bk^Ak 



2n (o + LOk) 



\UtBOF)+Y. 2 
^. (f2 + UJk) 



^^*Ak^*Bk 1+ .^.„ 2 



B ^kF) 



E 

k,l^k 



^*Ak^*Bl + ^*Al^*Bk 1^ ^ 11 \ 



(J] + Wfc)(217 + CJfc-|-W/) 



(22) 



is the result of dressing caused by interactions with both A and B together. 

In general, odd orders of the perturbative series have states with total spin (in 
the z direction) zero, and an odd number of phonons, whereas even orders have ±1 
total spin and an even number of phonons. 

2.3.2. Excitation of A and Selection of the Initial Dressed State The bare processes 



started with the initial state 



i'i ) = I^aIb Of) =(J^ |G(°)). What is the equivalent 

dressed state? In order to do that, let us operate with cr^ on |G). 

This is implemented using the approximation of an impulsive interaction, H' — 
aS (t) a^. Therefore, the evolution of the system neglecting the "free" parts is 

\^ {t = 0+)) = e-'""^ \G) (23) 
choosing a ~ 7r/2, the state is found to be (up to an irrelevant global phase) 

|V'(i = 0+)) = af |G) (24) 
hence we get that the "initial state", i.e. the state at t = is 



(1) 



(2,1) 



(2,2) 



(25) 



For simplicity, it is assumed here that Qa = = 
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where the corrections to the bare initial state 



V'f') = ITaIb Of) are 



(1) 



^*Ak HaIb Ifcf) + ^*Bk ItAta Ifcf) 



(26) 



and 



24)\ V^ (^Afc) + (^Bfc) 1^ I o \ I \ " ^Ak^Al + ^Bk^Bl |^ I 1 1 \ 

' ^ V2(i2 + Wfc)a;fc ^ 



kd^k 



{n + ujk) (wfe +0;;) 



(27) 



,,(2,2)\ \^ A^fcA_Bfc + X*Bk^Ak I I . „ V^A^fcA^fc I . o \ , 

/ = E 20 (1^ + .,) 1^^^^ + .kf ^ 

E ro!^'u9ot^'^l' J UTbUM (28) 

what can immediately seen from this result, is that the desired final state |4,At-B ^f) 
[jl is already contained in the dressed initial superposition. 

The next question is what is the correct approximation of the dressed initial state. 
There are several ways to answer this question: 

(i) Using |-0i) — (Jx \G) as an approximation for the dressed initial state and starting 
the time evolution. As mentioned before, in that case the probability we seek 
will not start from zero, since the final state already appears in the initial 
superposition. 

(ii) Measuring A's spin straightly after applying ^ and selecting only the up spins 
- which is like operating with ct^ instead; This should apply if one claims that 



Wi ) = 0"+ IG^"^) rather than V4 ) = 1^^"^) and therefore the same should 
be done with the dressed state. In this case, the effects of mutual dressing are 
lost (in the leading order) and hence the final state does not appear in the initial 
superposition and the probability starts from zero. This is the scheme used in [9J. 

(iii) Returning to the bare state: Measuring both A and i?'s spins, projecting to 
the bare initial state Ifyiis ^f)- This, of course, "ruins" the dressing process - 
however it is good for comparison. 

Finally, we note that in the dressed cases (the first two) the dressed states obtained 
perturbatively are not normalized. In order to normalize them (to order e^) one should 



add 



2 

^(0) \ ^ xhis additional term does not contribute to the calculations 



2 

and processes discussed in the following, and hence we ignore it 



2.3.3. Time Evolution of the Dressed State After selecting the initial state, it evolves 
until t — T similarly to the bare state. The time evolution of the dressed state is (we 

II We measure the bare states, since the measurement is assumed to take a very short time and is 
therefore unable to distinguish between the bare and dressed energies. 
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dt'Vit') t dt"Vit")\i,f^)~i f dt'Vit')\4''^) + \ipPyj +0{e^) 



The effect's amplitude is then (in leading order) 

^ (^) = ^'E i^*Bk^AkFik (t) + X*Ak^BkF2k (t)) 



where 



rfi'/o {t') e-^(f^+"'=)*' + 



2n {n + ujk) 



(29) 



(30) 



(31) 



F2k it) = - / dt'/o (f) e^t"-"")*' / dt'7o (i") g-*'^^-"^)*" 



dt7o (t') e*(f^-'^'')*' + 



idi 



^ + UJk Jo 



di 



2n{n + ujk) 



(32) 



and the constants di,d2 set the dressing (or non-dressing) scheme, as defined in the 
previous section: scheme 1 (cr;^) corresponds to di = 1, (i2 = 0, scheme 2 (u^) to 
di = 0,d2 = 1, and scheme 3 (no dressing) to di = d2 = 0. 

3. Examples 

3.1. Harmonic Chain 

As a first example, consider a linear harmonic chain with nearest-neighbour 
interactions. 

3.1.1. Properties 0} the Harmonic Chain Consider a circular (periodic) linear chain 
of A'" masses M, coupled by springs whose constant is described by the Hamiltonian 



n=l ^ 

As in |16j . we perform the canonical transformation 



(33) 



g„ = \ Mv\ l 



2K 



(34) 
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TT 



Pn = I = (35) 



and set 



2K 



2K 



(37) 



N-l 

hit) 



resulting in the Hamiltonian 

Hp = {pi + ql - ag„g„+i) (38) 

n=0 

Equation [37] implies 0<a<l. a— >lis the strong coupling limit, where 
— )■ oo. 

Periodic boundary conditions are imposed, and thus the equations of motion, 

qn + Elqn = Eq-^ (.qn+l + Q'ri-l) (39) 



yield (using £'0 — ) the dispersion relation 



^ n — /I — Q;cos6'fc 

t^fc = £^oV 1 - "cosfcife uW ; (40) 

V 1 — a 

and the solutions 

g„ (i) = ^ ^= fafce'^^'""-'"''*) + 4e-*(^'="-'^''*)) (41) 



where 



. iV — i I 

p„ (t) = _L ^ (a^e^^^'-"-"''*) - 4e-^(«''"-"''*)) (42) 

27rfc , , 



3.1.2. The continuum limit Let us take the limit N — >■ co, keeping the chain's 
length L constant: L = Nl — const. 
From the definition of a, one gets 

1 Ml/2 

- ^ 1 + ^ (44) 

Dimensional analysis shows us that the ratio M/K must be in units of [Time]'^. 
Therefore there exists a constant c, with velocity units, which satisfies 



(45) 
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it will be shown that this is the typical propagation velocity of the system. When 
taking the continuum limit, we shall also keep c constant. 
Now we get 

for small enough / (Z ^ 1, or equivalently L ^ N). 

(47) 

That gives the relation between N and a, which is dependent on the velocity c: 

" = (48) 

Therefore, for fc < A^, 



^fe« T\/47r2p + ^ (49) 



one can take y in'^k'^ + ~ -^Trfc for i/L ^ c (or smaller) 0. 
Hence, for small positive /c's, 

Wfc « 2TTky (50) 

3.1.3. Interaction with the spins and causality Denote the two oscillators which carry 
the internal spin as A and B (taking A and B as the oscillators' indexes), and define 
R = B — A as the angular distance, is units of Working again in the interaction 
picture, one gets the functions defined in (|14ll5p : 

F.(r).l;^£^^ii^^^^I^M) (51) 

k 



^-(-)"^E "°^"'"'''''^ (52) 

k 

By definition, uj^ — tON-k^ ai^d therefore 



■ N-l 



NFc{t,R) = — sin(wor)- 

UJQ 



tl 

V-2 



(53) 



^ Define the function g (k) = 2nk _^ jj^^ ^ _^ ]^ (.g^n be seen that for k = 1, the smallest 

yl+47r2fc2 k >oo 

interesting value of k (it will be shown why fc = is not of interest), g (1) Ki 0.987, so the convergence 
to 1 is very fast and the approximation is good even for small Ai's, as wanted. 
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and similarly 

^ cos(a)fcr-6ifcfl)+cos(aifcT+6ifcfl) odd N 



NF^ (r, R) = — cos {ujqt)+< 

UJo 



N -2 

A cos{uJkT-ekR)+coa(uJkT+6kR) , cos(u>n/2T-0n/2R) ^^^^ ^ 
^ /c=l 



(54) 

Let us analyze the time behaviour of the functions for a fixed R; The sums consist 
of trigonometric functions with coprimc frequencies, and tlierefore mostly do not add 
to some major contribution. However, there are special times in which that does 
happen. 

It is well known that the extremal points and zeros of sin u and cos u arc where 
u = for n € Z. Here it happens in two cases: UkT — OkR = ^ and LOkT+OkR = 
For the first case, the "special times" are 

= ^ + ^ (55) 
where the phase velocity of each mode is defined as 

^ft = -S- (56) 

using the small k approximation, 

therefore, for small k's, the phase velocity is not mode-dependent, although the 
frequency is. Therefore these "special times" are very near (almost equal) for small 
k's and n — 0. 

Looking at the zeros of each argument - times at which w^r — Oi-R = Oi one gets 
for small positive /c's: 

Tfe = — « — = - (58) 
Vk cN c 

where x = IR is the "real" distance along the chain, and hence c can be identified as 
the typical (and maximal) propagation speed along the chain, c: 

(59) 

The result is that several sine functions are equal to zero in very close times. 
Since in these zeros the arguments equal zero as well (turns from negative to positive, 
n = 0) - each of these sines is increasing, and hence they add up to a significant rise in 
Fc, around what is expected to be the "causal" time. This rise is not sharp, since the 
sum is not infinite; But it is significant. At higher fc's, for which the approximation 
does not apply, the arguments are zero at later times, but their contribution is less 
significant, because the zero points are less dense there (in time) and the functions 
are "weighted" with a ^ factor, getting smaller for larger values of k. 

The second significant rise occurs when uj^t + OkR = 0, again, in almost the 
same time for the small k modes. This time corresponds to the N — R distance (since 
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sin(ajfcT + 9kR) — sin(wfcr — 6^. {N — R))), which, remembering the periodic nature 
of the chain, represents the propagation in the other direction. 

What happens to Fa (r, R) at these times? When the cosine's argument is zero, 
it gets its maximal value, and it is symmetric in both the sides of the zero point: 
therefore, we can conclude that Fa (t, R) is almost symmetric for short times before 
and after the "causal" time! Therefore, the rise in the commutator "means nothing" 
for the anticommutator. 

Moreover, 

N-l ^ 

FjO,i?) =0- V — sin (6lfci?) = (60) 
k=i '^'^ 

(since 1/uJk is symmetric around k ~ N/2 and sin {9kR) is antisymmetric), as expected. 
So if Fc (0, i?) = and there is no significant rise (or fall) until r ~ ^, it would be 
reasonable that for earlier times the commutator is very small. 

This is analogous to a lightcone, but This does not imply exact causality, as in 
the continuum, since the commutator is not proportional to a step function 9 (t — 
However, it does show us that the commutator's contribution is insignificant comparing 
to the anticommutator's, for times "outside the lightcone", and hence the effect does 
not occur because of the nonvanishing commutator, or, in other words, the discrete 
nature of the system does not change the effect. An example for the behaviour of 
these "causality functions" is given in Figure [TJ It can be well seen from the graphs 
that the "lightcone" gets sharper as N grows - for smaller A^'s, the rise time is longer 
and hence the times for which the commutator is negligible are shorter than This 
is no surprise, and is due to the fact that as N grows the continuum limit, which is 
fully causal, is better approximated. 

Before turning to the example calculations, one should note that similar results 
should occur for discrete chains with other interactions, as long as the systems have 
translational invariance (and therefore the solution (|4ip is applicable, and one gets 
the functions (|51|) , ()52p ) , and their dispersion relation is approximately linear for small 




0.25 0.30 0.35 0.40 0.45 100 200 300 400 500 



Figure 1. On the left - Generation of a 'lightcone' - graphs of NFc (r, R = 0.3N) 
(as defined in Il52|l 1. for L = l,j/ = l,c=l;on the right - the commutator as a 
function of the distance R for A'^ = 1000, r = 0.31, L = 1, u = 1, c = I. 

3.1.4- Example Calculations - Bare States For example, fix the parameters N = 
100,L^l,nA = nB = 2,c=l,v = 1. 

The related commutator and anticommutator functions are plotted in figure [51 
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Figure 2. The magnitude of tlie commutator and tlie anticommutator as 
functions of the time, as defined in 1511521 1. for N = 100, R = 31, L = l,u = 
l,c = 1 

As a first example, we choose /a (t) = /s (t) = 1. Although this case forces the 
use of dressed states and does not allow the post selection of all three systems, we 
show these naive results as if these things were possible, in order to get the general 
(mathematical) picture. 

The results presented in Figure [3] are for two oscillators with a separation of 
i? = 31. 
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Figure 3. The probabiUty (on the left) and the magnitude of the separate 
contributions to the amplitude (on the right) of the transition |tA4-fl Op) — ^ 
14-AtB Of) as a function of time, for = 100, L = 1,Ua = = 2, c = 1, u = 

iJa (t) = Ib (t) = = 31. 

Next we shall show another example (Figure S]), the "physical" one, with the 
same parameters, except for the interaction functions, which are /a {t) = fs (t) = 
siii^ {jrf) 9 {t) 9 {T — t) : now the interactions are limited in time, so if we take T short 
enough, dressing effects do not have to be considered - and therefore this is a more 
exact result than the previous example. In the example shown in the following figure, 
T = 0.1 which is smaller than the "causal time" T ~ 0.31. Moreover, turning the 
interactions off allows us to perform the "nonlocal measurement" afterwards. The 
results are shown in figure H) 
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0.04 0.06 



Figure 4. Graphic Results, showing the total probability of the transition 
ItA^B Oi?) — > 14-Ats Op), for the case of time-limited interactions: N = 100, L = 
1,0^ = Ofl = 2,c = 1,1/ = 1,/aW = /flW = sin^^) e{t) 9 {T-t),T = 
0.1, -R = 31. In this case, the commutator's contribution Ac is highly negligible 
comparing to Aq. 

3.1.5. Example Calculations - Dressed States In order to get the dressing effects, the 
interaction functions were chosen to be 

Ia (t) ^ fs (t) ^ \ 2 ( ■nt \ a u\ a Irr j.\ ^ \ r, ^^^^ 

[/o(t) =cos^(^)6'(t)6l(T-t) t>0 

for a very large r. It corresponds to starting with a dressed ground state at t = 0, 
and tlren switching the interaction off after a duration T. 

Plugging the harmonic chain's eigenvectors into ([28|) . we get that the time 
independent part of the amplitude for harmonic chains is 

'^^^^ mnujk (n + cjk) ^^^^ 

k 

defining, as usual, the angular distance R = B — 

V^V COs(gfci?) 

this is, of course, a dressing effect. If the function G {R) for a given chain with N 
oscillators is examined, it can be seen that it is a decreasing function: the "zero-time" 
amplitude is bigger for closer oscillators. The minimal value, however, is not zero: 
dressing occurs for every distance R. 

The maximal distance is reached when the oscillators are on antipodal sites on 
the ring: i.e., Rmax — -y, for a chain of an even N. Therefore one can define 

= i:«TTiy 

noting that 

2TrkN 

OkRmax = = (65) 
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one gets 

(for even N). As one can see in Figure [Sj this function decreases rapidly as the chain 
gets longer. 




100 200 300 400 500 2000 4000 6000 8000 10000 



Figure 5. -i^G as a function of R, for = 1000 (on the left), and -^Gmin (N) 
(on the right), for even values of N. L = 1,Q = 2, c = 1, v = 1. 

Using p0|31l32p . the total amplitude is given by 

^ (*) = ^T. - {^"'""'''F^k + ^^''"''■^2. {t)) (67) 
iiy — Ulk 

k 

We show two example calculations - as before, the first one is for general 
impression and has /o {t) = 1 (Figure [6]) and the second one, which is more relevant, 
has /o {t) = cos^ (^) e (t) 0{T-t) (FigureEl). In both the examples, N = 100, L = 
1, r2 = 2, c = 1, = 1,R = 31, and hence the "causal time" is T = 0.31. In the second 
case, T ~ 0.1, which is smaller than the causal time, as expected. One can see that in 
the "real" interaction, where the interactions are switched off (Figure[7]), the measured 
probabilities (i.e. in the end of the interaction time) of the dressed processes are much 
larger than in the bare ones; Moreover, when the first dressing scheme (with a^) is 
used, the probability is higher, in our example, in about two orders of magnitude than 
in the second scheme (<t+), and is much more than in the bare case. For such a short 
time the constant contribution to the probability, resulting from the appearance of the 
desired final state in the initial dressed superposition, dominates, and hence, as can 
be seen in the figure, the probability of the first dressing scheme is almost constant in 
time. 
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Figure 6. Graphic Results, showing the total transition probability in the three 
dressing cases, for N = 100, L = l,n = 2, c = 1, u = 1, R = 31, fo (t) = 1. (Pi - 
the initial dressed state is generated by operating with on the dressed ground 
state; P2 - with cr^; P3 - Bare initial state, no dressing at all.) 
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Figure 7. Graphic Results, showing the total transition probability in the three 
dressing cases, for N = 100, L = 1,Q = 2,c = l,u = 1,R = 31, /o (t) = 
cos^ (^) e (t) e (T - t) ,T = 0.1. (Pi - the initial dressed state is generated by 
operating with on the dressed ground state; P2 - with a^; P3 - Bare initial 
state, no dressing at all.) It can be seen that the probability in the first dressing 
scheme is much larger than the other two; This is no surprise since the desired 
final state is included in the initial dressed superposition. As mentioned, the 
"causal time" is T = 0.31, and it can be seen that since this dressing (constant) 
contribution dominates the amplitude before that time, the probability is almost 
constant. 

3.2. Ion Trap 

In this section we study the Fermi Problem for trapped ions. 

Consider a linear ion trap, consisting of N identical ions in a Linear Paul trap. 
The ions' free linearized Hamiltonian, Hp amounts to a set of free harmonic oscillators 
(as in (12)); The displacement (coordinate) operator of each ion is 
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where a;fe,afe,aj. are as usual, and D„fc, elements of the eigenmode vectors, are 
calculated numerically (as well as the frequencies w^). In the Lamb-Dicke limit, using 
rotating wave approximation for the internal degrees of freedom only and moving to 
the rotating frame, we have 

H' = -^d{a^ + af)+HF + e ^ /„(i)^:g« (69) 

ne{A,B} 

where 6 = uJi^ — Q is the detuning from the laser frequency, thus one gets, in interaction 
picture 

HL^^ E /«(t)Ke-^**+a!!e^^*)<z„(i) (70) 

ne{A,B} 

which is similar to the general interaction Hamiltonian presented before for discrete 
systems: one has only to replace 17 with —5. 

The perturbative calculation turns out to be very similar to the harmonic chain 
case. We therefore proceed with a nonperturbative example. 



3.2.1. A Non Perturbative Example As seen before, although the effect exists, it is 
very small. In order to increase its magnitude, one would like to increase the coupling 
constant e; and use instead a non perturbative approach. The effect can then be 
anticipated to become significantly larger. 

We demonstrate this by following and considering as a "proof of principle" 
an ion trap with only two ions. Their ground state, \V) = |00) (in the eigenmodes 
basis), can be written in the local number basis as 

oo 

\V)^Zj2e-P-\nA)\nB) (71) 

where 

Z = (72) 

It is a thermal state, whose "temperature" /3 can be calculated, for example, using 
symplectic diagonalization and Williamson modes fW] . 




where A is the symplectic eigenvalue 



and ujq^uji are the frequencies of the system's eigenmodes. 

Calculating the frequencies of our system (wi = \/3ujo), we get A « 0.5189, 
and hence /3 w 1.9916, ~ 0.1364. Due to the smallness of e^^, one can write, 
approximately, 

\V) « \v) = \0a) \0b) + \1a) \Ib) (75) 
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The displacement operators, in terms of the local sites creation and annihilation 
operators, are given by 

Qn = (On + 4) (76) 

T, the time during which the interactions are open (the pulse length) must satisfy 
T ^ ^ - shorter than the typical propagation time along the trap, and thus it obeys 
the "outside the lightcone" demand: for ion traps, 

P / sr--DakDbk cos {ukT) 

^« W = (77) 

and 

F,(r) = ^''"^-'^"'""("^") (78) 
k '^'^ 

and so, for T < T <C one gets that (r) is negligible comparing to Fa (r). 

We consider the case that the interaction coupling constant is sufficiently strong, 
using strong and impulsive laser pulses, such that the free evolution Hamiltonian can 
be neglected. One can then approximate, H Hi, and get 

\iP {t)) « e<"^'^^("+"*)+"«'^» IIaIb Vf) (79) 

where a, act on A and 6, fe^ on B, and q:„ = — -^== dtfn (t) . 

Working in the Schrodinger picture, the operators are time-independent. 
Therefore the time evolution operator can be separated into two commuting parts 

1^ {t)) ^ e*"^"^("+"')e^"«"-''(''+''') ItAiB Vf) (80) 

which can be decomposed further: 

and so one obtains 

IV' (i)) « e-^^^e^«^'^^«'e^«^""'''e'«^'^^^«e^"^'^"'' I^aIb) (|0a) \0b) + |U) |1b)) 

(82) 

note that 

^iaAa^a^ias.Sb |^^^^) |0^) |0s) = jtA^s) \0a) \0b) (83) 

and thus 

e""^ (1 + iaA(7^a'< + •••)(! + iaea^b^ + ■■■) ItAis) \0a} \0b) = 

- e-^^QAas lUts) I 1a) |1s) + ■4') (84) 



Similarly, 
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|tAiB)|U)|lB) = 

- e ^^aAae lUts) |0a) IOb) + ip' 
where (jpfl'if^ = (^ipfl'ip''^ = 0. The AmpUtude is 

A (q;a, as) ~ — 2e = ay^ase"'^ (85) 

and therefore 

P {aA,aB) ~ 4e-"i+"«a^a|e^2/3 (gg) 

if the similar pulses are used for both the ions (a a — o-b = ck), the transition 
probability is 

P(a)«4e-2"'aV2'3 (87) 
This probability is maximal for a = 1: P (1) w e-2/3-2 ^ 0.0100819 w 1% . 

4. Summary and Discussion 

We have discussed the simulations of the Fermi two-atom problem within the context 
of discrete systems. The fact that in such cases the interactions between the "atoms" 
and the "field" are controllable, enables us to perform post selection on all the 
participating systems, and therefore the problem could be re-examined in Fermi's 
original formulation [T]. 

As shown in section 2, the general calculation of the relevant transition amplitude 
for a discrete system, to leading order in perturbation theory, can be decomposed into 
two separate parts - one is correlation-dependent, similar to the continuous case outside 
the lightcone, and the other is commutator-dependent. The second part vanishes in 
the continuous case, and hence defines a light-cone analogy for our system: times 
at which the second part is negligible comparing to the first are defined as being 
"outside the lightcone", and then the first part dominate and the discrete effect is 
similar to the continuous one. We have seen, for example, for harmonic chains, that 
the commutator's contribution grows around t ^ ^, and as the number of oscillators 
grow, the discrete "lightcone" is approaching a real, continuous lightcone. 

Another thing of interest was the question of dressing. Turning the interaction 
on quickly and for a short time enables us to use the bare states as the physical 
states of the system. However, if the interaction is turned on adiabatically, the states 
undergo a dressing process, and the physical states are the dressed ones. The dressed 
states were obtained perturbatively (as in |151 19]). Two different dressing schemes 
were considered, and in one of them the desired final states already appears in the 
initial dressed superposition. Moreover, even in the bare case the nonlocal excitation 
cloud arises immediately, as discussed in appendix 2. 

The effect's probability is very small. In order to obtain measurable probabilities 
one could try to increase the interaction strength, leaving the perturbative regime. We 
have shown, as a proof of principle, an example for an ion trap of two ions with a very 
strong, impulsive interaction, for which a probability of about 1% was calculated. 
This, however, was done only with the bare states, since our dressing process was 
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perturbative and hence dependent on the interactions being weak. In an experiment, 
however, one could turn on strong interactions adiabatically, causing the system to 
dress, and observe the influence on the measured probability. 

The field's initial state could be considered further as well: for large N, it might 
not be possible in current day experimental technologies to cool the system to its 
ground state. Hence the effect with the field in a thermal state of a low temperature 
is another question worth addressing. 

In conclusion, one can see, in a discrete system for which a lightcone can be 
approximated, that the probability of swapping the atoms' internal degrees of freedom 
appears to be noncausal as in the continuous case. As in the continuous case, there 
is a nonlocal "cloud" of field excitations which probably cause this to happen. The 
discrete effects, with the benefits which result from turning on and off the interactions, 
can serve as a quantum simulator of Fermi Problem, showing the effects of dressing. 
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Appendix 1: Perturbative Calculation of the Bare Amplitude 

Working in the interaction picture, we get (assuming that e is small enough) that at 
time < > 0, 

1^1 - I jj* dt'Hi {t') - dt'Hi (t') dt"Hi [t") + O (e^)^ \i> (0)) (A.l) 

Considering the relevant (bare) initial and final states (|6l7p . the amplitude we 
seek is: 



O {e') \tAiB Of 



A (t) = {UtB OfI ( 1 - i / dt'Hi (f) - I dt'Hi (*') / dt"Hi {t" 
\ Jo Jo Jo 

(A.2) 

The lowest order process contributing to this amplitude is an exchange of a virtual 
phonon between the two "spins". This can be written as 



t' 

3\ 



A (t) - - {UtB Of I / dt'Hi (f) / dt"Hi it") ItAis Of) + O (e^) (A.3) 
Jo Jo 

Plugging in the interaction Hamiltonians ([5]) the amplitude can be decomposed into the 
two possible orders of this process (emission - absorption and absorption - emission) 

A{t) = - {UtB Of I /* dt'Hf it') f dt"Hf {t")+ f dt'Hf {t') f dt"Hf {t") {^aIb 0f)+O (e') 
Jo Jo Jo Jo 

(A.4) 



After operating on the spins' spaces, the amplitude (in lowest order) is 



The Fermi Problem in Discrete Systems 



22 



Ait) 



t i^r 



^2 = j^dt'fA{t')e-'''-' dt"fB{t")e'''-' {G\qA(t')qB{t")\Q) + 
dt"fs{t")<^^-'" f dt'fA it') e-"-^^' (0| qn it") QA (f) |0) ^ - "^^^ + 



^0 



(A.5) 



Changing the order of integration and the order of operators in both the terms 
Aab it) , Aba it) we get 



Aab it) 



= f dt"fB it") e'"-*" /* dt'fA it') e-'"-*' (0| qB it") qA it') |0) 
Jo Jt" 



+ f dt"fB it") e'"-*" /* dt'fA it') e-^"-^*' (0| [qa it') , QB it")] |0) (A.6) 

^0 Jt" 



^ /* dt'fA it') e-^"-**' /* dt"fB it") (0| qA it') qB it") |0) 



'0 Jt' 

dt'fA it') e-'"-*' / dt"fB it") e'"«*" (0| [qB it") , qA it')] |0) (A.7) 
Jt' 

adding each of the changed terms to the other term unchanged, a new expression for 
the amphtude is reached. Adding the two expressions obtained that way and dividing 
the sum by 2, one gets 

Appendix 2: Estimation of the Field Excitations Distribution 

Another thing of interest is the distributions of field excitations ("phonons") in the 
system, as a function of the sites. Following ^9J, we define the excitations distribution 
at site n, time as 

A^W = (^.(<)|g+(<)g-WIV^.(t)) (B.l) 

(in interaction picture), where 

N-l 

<lt{t)=Y.^*nAe'^'' (B.2) 

N-l 

q- it) = Y^Kkake-''^"' (B.3) 

create/ destroy a field excitation at site n. 

Dn it) - ^A;,.A„,e^("^-"')* it)\alai |V^, it)) (B.4) 

k,l 

using (Uni), 

A. it) = e2^A;fcA„,e*(-'=-')* ic*A, it) cm it) + c*b, it) cbi it)) + O (e^) (B.5) 

k,l 
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where 

CAk it) = (^Y^ + * 1^ dt'fo {t') e-(^-'=)*') (B.6) 

CBk (t) = X*B, + ^ 1^ dt'fo it') e^("+-'=)*') (B.7) 

di , d2 are the dressing factors defined previously. It can be seen, therefore, that it the 
initial state was bare, D„ (t = 0) = as expected - but the "tails" of field excitations 
will be spread out for any t > 0\ 

As before, the leading order does not involve any mutual dressing factors. In 
other words, in the leading order, the phonon distribution is just an addition of the 
field excitation distributions of the two sites coupled to the field. Therefore we can 
deduce what is the field excitations distribution of a single site: if it is in a spin-up 
state (such as A), 

Dl (t) = e^^A:,A„,e'(-'=-"')*c^, {t) cm (*) + O (e^) (B.8) 

k,l 

and if it is in a spin-down state (such as B), 

Di (t) = e2^A;,A„,e*('^'=-'^')*cj3, {t) cbi (t) + O (e^) (B.9) 

k,l 

and again, for the Fermi Problem, in leading order, 

D^{t) = Dl{t)+Di{t) + 0{e'') (B.IO) 

One should note that despite the fact that q^Qn could serve as a qualitative 
analogy to the local site phonon number operator, it is a nonlocal observable since it 
cannot be expressed by the local operators qn,Pn- Therefore it would not be possible 
to simulate the results of this appendix, and it could serve only for a qualitative 
impression. 

Example for the Harmonic Chain We present an example for the generation of 
"nonlocal cloud", when the initial state is bare (figure Bl). 
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Figure Bl. A Plot of the local field excitations distribution, for an initial bare 
state, of the single spin-up atom in the 50th site of a chain of 100 sites, with A'^ = 
100, L = 1, O = 2, c = 1, 1/ = IJa it) = Ib (t) = sin2 [^) e (i) e {T - t),T = 0.1. 
The generation and spreading of the nonlocal "cloud" can be seen. 
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